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Abstract. Gamma-ray bursts (GRBs) must originate from low baryon 
load, ultrarelativistic outflows; however, slower, more baryon-rich out- 
flows (BROs) should also arise in connection with GRBs as "circum-jet 
winds" and/or "failed GRBs" . We study the possibility of nucleosynthesis 
within such BROs by conducting detailed reaction network calculations in 
the framework of the fireball model, showing that they can be interesting 
sites for synthesis of heavy neutron capture elements, as well as of light 
elements such as deuterium. These products may be observable in the 
companion stars of black hole binary systems or in extremely metal-poor 
' stars, offering an interesting probe of conditions in the central engine. 
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1. Introduction 



Successful generation of GRBs requires the formation of an ultrarelativistic out- 
flow with bulk Lorentz factor V ^ 100, implying very low baryon-loading (e.g. 
Meszaros 2002). Since the temperature at the base of the outflow should be 
of order MeV and the baryons are likely to contain a high fraction of neutrons 
(Pruet, Woosley & Hoffman 2003, Beloborodov 2003), some production of nuclei 
starting from free protons and neutrons is expected to occur in the expanding 
flow. However, as shown in several recent papers (Lemoine 2002, Pruet, Guiles 
& Fuller 2002, Beloborodov 2003), in the very high entropy and extremely rapid 
expansion characteristic of GRB outflows, nucleosynthesis is limited to only 
small production of D and 4 He, which are difficult to observe. 

Nevertheless, most GRB progenitors should also give rise to associated out- 
flows with higher baryon- loads and lower velocities (baryon-rich outflows, or 
BROs). It is highly probable that different types of "circum-jet winds" of baryon- 
loaded material surround and coexist with the narrowly collimated GRB jet. 
The GRB jet production mechanism is likely to act not only on the baryon-poor 
zones near the jet axis, but also on the more baryon-contaminated, peripheral 
regions, leading to BROs as a natural byproduct. Models involving core collapse 
and jet penetration in massive stars (MacFadyen and Woosley 1999, Wheeler et 
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al. 2002) can also induce BROs through entrainment and mixing with the stel- 
lar material. In black hole accretion disk models, baryon-rich winds may arise 
from the outer disk by mechanisms distinct from the GRB jet (Janka & Ruffert 
2001, Narayan, Piran & Kumar 2001, MacFadyen 2003). Some observational 
evidence support the existence of such circum-jet BROs (e.g. Woosley, Zhang 
& Heger 2002). Alternatively, BROs may occur without concomitant GRBs 
as "failed GRBs" , potentially with event rates higher than successful GRBs, if 
the baryon-loading process acts more thoroughly or the the GRB jet-driving 
mechanism operates less efficiently in the central engine. These may possibly 
be identified with observed "hypernovae" (Nomoto et al. 2002), X-ray flashes 
(Kippen et al. 2002), or with some hitherto unrecognized type of transient. Ei- 
ther way, the lower entropy, slower expansion and higher ejecta mass of BROs 
compared to GRB jets make them much more interesting from a nucleosynthesis 
viewpoint. As a first study of nucleosynthesis in BROs associated with GRBs, 
we investigate this problem utilizing the simplified dynamical framework of the 
basic fireball model, but incorporating detailed nuclear reaction networks in- 
cluding both light and heavy elements. More details can be found in Inoue et 
al. (2003). 



2. Nucleosynthesis in GRB-BROs 

We consider a wind fireball, i.e. an adiabatic, thermally-driven and freely- 
expanding outflow (Paczyhski 1990). The key parameters are the luminosity 
L = 10 52 Z/52erg s _1 , the initial radius r$ = 10 7 rojcm, and the dimensionless en- 
tropy rj = L/Mc 2 . The fireball first undergoes an acceleration phase, where the 
comoving temperature and density decrease exponentially on the initial dynami- 
cal timescale too — 0.33 msec roj. After the internal energy is converted to bulk 
kinetic energy, there follows a coasting phase where T and p b decrease as power- 
laws, asymptotically becoming T = To?7 _1 (t/t ( io) ^ 3 and pj, = pb,o r l 3 (t/tdo) ' 2 - 
These formulae are valid for rj ^> 1. Although r] <^ 1 (nonrelativistic expansion) 
may be realized in some progenitor models, here we assume rj > 2, corresponding 
to baryon load mass M b < 2.8 x 1O~ 2 M £^5 3 for total energy E = 10 53 erg£ , 53 , 
which may be a reasonable upper limit for relatively baryon-free progenitors 
such as neutron star mergers (Janka & Ruffert 2001) or supranovae (Vietri & 
Stella 1998). The initial electron fraction Y e is taken to be an additional pa- 
rameter, in view of the uncertainties in its expected value (Pruet et al. 2003). 
We implement nuclear reaction network codes including a large number of light, 
very neutron-rich nuclei (Terasawa et al. 2001). 

Taking fiducial values of L52 = 1 and tqj = 1, we concentrate on the effect 
of variations in rj and Y e . Although 4 He is the dominant final species besides 
remaining free nucleons for all cases, interesting amounts of other light elements 
can also result, particularly deuterium, whose final mass fraction can reach Xd ~ 
0.02. (Fig.l). For r] ;> 20, the reactions mainly occur during the exponential 
expansion (acceleration) phase, similar to the case of big bang nucleosynthesis. 
Contrastingly, for i] <; 10, nucleosynthesis proceeds mainly during the power- 
law expansion (coasting) phase. When Y e <^ 0.5, free neutrons remain abundant 
in the outflow until they start decaying into protons, so that D production can 
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Figure 1. Contours of final abundance by mass fraction of deuterium 
for different parameter values of n and Y e . 

proceed by subsequent p(n,7)d reactions at very late times. Other elements such 
as Li, 9 Be and n B can be generated at levels of X ~ 1CP 8 — 1CT 7 . 

The synthesis of heavy elements through neutron capture is even more in- 
triguing. The fiducial n = 2 outflow is characterized by very high entropy 
(s/ks — 2500), low density (pb,o — 1-5 x 10 5 g cm~ 3 ) and short initial expansion 
timescale (i^o — 0.33 ms), leading to nucleosynthesis which is markedly differ- 
ent from the better-studied r-process in supernova (SN) neutrino-driven winds 
(e.g. Qian 2003). The final abundances are shown in Fig. 2. When Y e <J 0.4, 
considerable reaction flows can occur to high mass numbers, although the peak 
abundances are at levels of Y ~ 10~ 6 . (Note that the apparent peak at A ~ 250 
is an artifact of our neglect of fission.) The abundance patterns do not match 
the solar distribution and possess peaks in between the solar s- and r-process 
peaks. This reflects the occurrence of reaction flow paths intermediate between 
those of the s- and r-processes, resulting from the balance of neutron captures 
and /3-decays rather than (n,7)-(7,n) equilibrium. Such processes have been dis- 
cussed in previous contexts and called the "n-process" (Blake &; Schramm 1976). 
We also note that the neutrons remain abundant after the reactions freeze out 
until they freely decay, instead of being exhausted by neutron captures. 



3. Observational Implications 

The present results of nucleosynthesis in GRB-BROs have interesting observa- 
tional implications. For our case of r] = 2 and Y e ~ 0.1 — 0.4, the ejected mass in 
heavy elements per single GRB-BRO event should be M^^RO ~ 1O _7 -1O _6 M0 
for A ^ 100, comparable to or somewhat less than that expected in r-process 
elements from a normal core collapse SN (Qian 2003). Considering the fact that 
the abundance patterns are different from solar, and that the intrinsic event 
rates of (either successful or failed) GRBs are estimated to be less than that 
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Figure 2. Final abundances as function of mass number, for r/ = 2 
and Y e = 0.1, 0.3, 0.4, 0.48 and 0.498, represented by the dotted, solid, 
dashed, long-dashed and dot-dashed curves, respectively. The upper- 
most, thin solid curve is the solar abundance distribution in arbitrary 
units. 

for SNe, the contribution of GRB-BROs to the solar /Galactic abundances are 
unlikely to be large, at least for the conditions treated here. 

However, some observable environments may manifest the local abundance 
pattern of the BRO ejecta. Many facts suggest that chemical evolution in the 
early Galaxy proceeded inhomogeneously, and also that GRBs are associated 
with the most massive progenitor stars. If so, these explosions may have dom- 
inated at the earliest epochs, so that the BRO nucleosynthetic signature are 
visible in the most metal-poor stars. Alternatively, both light and heavy ele- 
ment products may be clearly discernible in the companion star of a black hole 
binary system, where the event that formed the black hole may have left traces 
of the explosion ejecta on its surface (Israelian et al. 1999). High resolution 
spectroscopic observations of such stars are warranted for detailed studies of 
their elemental abundance patterns. Such studies should yield important con- 
straints on the value of Y e , providing us with valuable insight into the innermost 
regions of GRB sources which are otherwise very difficult to probe. 

4. Conclusions and Outlook 

In this exploratory study, we have shown that BROs associated with GRBs can 
be interesting sites for nucleosynthesis of both light and heavy elements, with 
important observational consequences. We have neglected here various realis- 
tic effects that can be potentially important, including fission and proton-rich 
nuclei, collimation and higher baryon loading, neutrino irradiation, nonthermal 
spallation, neutron captures on external matter, etc., all of which can be sig- 
nificant, e.g. in the widely discussed collapsar scenario. Until we arrive at a 
definitive understanding of the physics of the GRB central engine, construction 
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step by step of increasingly realistic models incorporating these effects and com- 
parison with various observations may help us constrain the unknown properties 
of the GRB source, as well as to clarify the role of GRB-BROs in the chemical 
evolution of the Galaxy It is certain that much interesting work lies ahead. 
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